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The sympathetic cooling of an initially degenerate Fermi gas by either an ideal Bose gas below T c 
or an ideal Boltzmann gas is investigated. It is shown that the efficiency of cooling by a Bose gas 
below T c is by no means reduced when its heat capacity becomes much less than that of the Fermi 
gas, where efficiency is measured by the decrease in the temperature of the Fermi gas per number of 
particles evaporated from the coolant. This contradicts the intuitive idea that an efficient coolant 
must have a large heat capacity. In contrast, for a Boltzmann gas a minimal value of the ratio of 
the heat capacities is indeed necessary to achieve T = and all of the particles must be evaporated. 

PACS numbers: 



I. INTRODUCTION 

Sympathetic cooling of fcrmionic atoms has proven to 
be an efficient tool to produce strongly degenerate Fermi 
gases pi 0, 0, 0, • It is often thought that when the 
specific heat of the coolant becomes less than that of the 
fermions the cooling process slows drastically and even 
comes to a halt 0, 0, ■ This belief derives from the fol- 
lowing argument based on standard thermodynamics 0] . 
A coolant is not only described by its temperature but 
also by its heat capacity. The equilibrium temperature 
between two bodies initially at different temperatures Tf 
and T b is T = (C f T f + C b T b )/(C f + C b ), where Cf and C b 
are the heat capacities of the two bodies 0. Therefore, 
regardless of having a very cold object b (T b <C Tf), it 
will not act as a good coolant unless it is also true that 
C b » C f . 

However, this argument does not include the evapora- 
tion of atoms from the coolant 0,^1, a key ingredient of 
sympathetic cooling. In the following, we shall consider 
the ideal case of sympathetic cooling, where all losses of 
particles due to inelastic collisions are absent. One can 
then evaporate sufficiently slowly that the Fermi gas and 
the coolant remain close to thermal equilibrium. This is 
in the spirit of the above argument, based on heat ca- 
pacities, as heat capacity is a property of an equilibrium 
system. In this case, the only limitation of resources 
is the number of particles in the coolant. We then in- 
troduce an idealized model of sympathetic cooling: the 
atomic motion is treated in the semi-classical approxi- 
mation and cooling proceeds in discrete steps. In each 
step, the atoms of the coolant go through an evaporation 
phase, and then the fermions and the coolant thcrmalizc 
(see Fig. n below). This is a generalization of the dis- 
crete scheme used for evaporative cooling in Ref. to 
the case of sympathetic cooling. A similar semi-classical 
and discrete model for sympathetic cooling was put for- 
ward in Ref. , with the addition of features like colli- 
sional losses and space dependent effects in the evapora- 
tion process, but without a determination of the limits of 
the cooling process. Using our idealized model, we shall 
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FIG. I: Shown is a diagram of the evaporative cooling process 
in our model. In the upper two panels, a gas of Nf fermions 
at temperature T is in thermal equilibrium with a coolant of 
Nb bosons or Boltzmann particles. In the middle two pan- 
els, the tail of the number distribution of the coolant is cut 
and SNb particles are removed. In the lower two panels, the 
systems are permitted to rethermalize, resulting in a reduced 
common temperature T — ST. The values of ST and SNb have 
been exaggerated for the purposes of illustration; the units 
are arbitrary. 



then show that the argument based on the heat capac- 
ities is approximately correct for fermions cooled by an 
ideal Boltzmann gas 0], but incorrect when cooled by 
an ideal Bose gas below the Bose condensation tempera- 
ture T c . In fact, in the latter case, it is preferential to 
have C b <C Cf to minimize the number of particles of 
the coolant whose evaporation is required to reach zero 
temperature [l5j . 



II. THE MODEL 

We assume that the Fermi gas, denoted by the sub- 
script /, and the coolant, denoted by the subscript b, are 
initially in thermal equilibrium at a common tempera- 



ture T. One evaporates all particles of the coolant with 
an energy greater than rjk B T 0,0,0,0, where r\ is 
a constant, typically 5 to 7 in current experiments. This 
has the effect of removing a number of particles SN b and 
an energy SE b from the coolant: 



5N b = / dep b (e)N b (e), 

Jt)k B T 

5E b = / dep b (e)N b {e)e, 

Jr)k B T 



(1) 

(2) 



where p b (e) is the density of states of a particle in b and 
N b (e) is the mean occupation number of a state of en- 
ergy e. In what follows we restrict ourselves to the semi- 
classical regime in which the number of trap levels with 
an energy less than k B T is much larger than unity, both 
for species / and b. Then the case of a density of states 
of the form 



PfA e ) = A f,b 



: S+l/2 



(3) 



may be considered: this encompasses both a harmonic 
potential (S = 3/2) and a box ((5 = 0) Here the 

subscript / (6) stands for the fermions (coolant). In 
the two limiting cases of a non-degenerate gas N b (e) = 
cxp[— [3{e — p)] and a Bose gas below T c . N b (e) = 
l/[exp(/3e) — 1], the mean energy of an evaporated par- 
ticle is independent of the chemical potential p and 
proportional to T: 



5E b 
SN b 



f(v)k B T. 



(4) 



The function f{rj) is different for a Bose gas below T c and 
for a Boltzmanngas, but has the same large rj behaviour 
in both cases [lg : 



/(r?)=r?+l + 



(5) 



Once the particles are removed, the system is left to 
rethcrmalizc. The decrease in temperature ST is ob- 
tained by conservation of energy 

E f {N f ,T)+E b (N b ,T)-SE b 
= Ef (Nf , T — ST) + E b (N b — 5N b , T — ST) . (6) 

Note that Ef is the total energy of the fermions, not 
the Fermi energy. The parameter r\ is chosen sufficiently 
large so that SN b < N b and ST «T. Equation © may 
be expanded to leading order in SN b and ST: 



dE b 
dN b 



(N b ,T) 



ST 



SN b C f {N f ,T) + C b (N b ,T) ' 



(7) 



where the heat capacities are defined by C/ = dxEf and 
C b = drE b . 

Equation Q, valid for a Boltzmann gas or a degenerate 
Bose gas, will be studied for both a three dimensional 




FIG. 2: Shown is the temperature of a degenerate Fermi 
gas sympathetically cooled by a Bose gas with a condensate 
present, as a function of the number of bosons evaporated in 
a harmonic trap, that is of the difference between the initial 
number of bosons Nb,o (assumed to fulfill Ea. l|14|l ') and the 
final number iV;,. Two cases are depicted: uif/uib — 1 (solid 
line) and ujf/ujb = 5 (dashed line) where ujf/uib is the ratio of 
the mean trapping frequencies of the fermions and bosons. It 
is apparent that Uf/uob ^5 1 requires a greatly reduced num- 
ber of bosons in order to cool to T = 0, and may therefore 
be considered as a more efficient experimental configuration. 
Here the choice of parameters was 77 = 6 for the energy cut-off 
and To/Tf = 0.2, where To is the initial temperature. 



box and a harmonic potential. In the former case, the 
mean field energy per particle of the coolant, due to the 
interaction with the Fermi gas, is constant and therefore 
plays no role. In the latter case, it is not exactly constant, 
but we shall assume that it is negligible as compared to 
k B T . Additionally for simplicity it will be assumed that 
the initial state of the fermions is degenerate. In this case 
the total energy is 




E f (N f ,0) 

41 



(8) 



where a = tt 2 (8 + 3/2)/6, k = tt 2 (5 + l/2)(<5 - 3)/20, 
S is the exponent in the density of states and the Fermi 
temperature Tp is given by 

A f (k B T F ) 5 + 3 / 2 = (5 + 3/2)JV>. (9) 

In Eq. JHJl we shall keep terms only to second order in T 
in the following. 



III. COOLING BY AN IDEAL BOSE GAS 
BELOW T c 

A non-interacting Bose gas with a temperature below 
T c has a total energy 

vA b 



E b (N b , T) 



5 + 5/2 



(k B T) 



5+5/2 



(10) 
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where v = (5 + 5/2)r(<5 + 5/2)C(5 + 5/2) is a constant and 
r and C are gamma and zeta functions, respectively |l9j |. 
Note that Eb does not depend on the total number of 
bosons, since the occupation of the excited states is sat- 
urated to its maximal value in the presence of a conden- 
sate. Therefore dEb/dNb = and Eq. reduces to 



6T 
5N~b 



f(v)T F 



2aN f 



1 + 



3v Ah 



TT 2 Af \T, 



T 



5+1/2 



(11) 



where we used the fact that Af is related to the Fermi 
temperature Tp by Eq. 0. When species / and b are 
in a box, A b /Af is simply the ratio Vb/Vf of the spatial 
volumes occupied by the two species. When species / and 
b are trapped harmonically, Ab/Af is given by the ratio 
(uif /ci>b) 3 where uif,b are the geometric means of the three 
oscillation frequencies for species / and b, respectively. 

The term in T s+1 / 2 in the denominator of Eq. ltTT|) 
represents the ratio of the two heat capacities: 



c f 



Zv A h 

IT 2 Af 



T 



5+1/2 



(12) 



If one measures the efficiency of the cooling by the tem- 
perature decrease per number of bosons evaporated, i.e., 
ST/5Nb, it is most efficient to have Cf, <C Cf. In this 
case the first term in the denominator dominates and the 
derivative is maximized. In Fig. [21 is shown the temper- 
ature as a function of the number of bosons evaporated, 
as can be obtained from the solution of the differential 
equation Eq. (|11|) for the case of harmonic trapping and 
r\ = 6. It is apparent that it is advantageous to evapora- 
tive cooling to choose a trap frequency ratio u>f/u>b ^ 1, 
in order to maintain Cb/Cf small at typical experimental 
temperatures T/T F ~ 0.2 or 0.1. 

A simple condition may be obtained from Eq. I|ll|) for 
the total initial number of bosons necessary to cool to 
T = 0: 



N b n > AN b = 



C 



f,0 



k B f(v) 



1 



1 



Cb.o 



5 + 3/2 C /i0 



(13) 



where T is the initial temperature, Cq are the initial heat 
capacities and AAf, is the number of bosons evaporated. 
As the coolant is assumed to be initially a Bose gas with 
a condensate present, the initial number of bosons must 
exceed the saturated value of the population of the ex- 
cited states: 

N b , > A h T(5 + 3/2) ((6 + 3/2) (k B T Q ) 5+a / 2 . (14) 

Note that, as an arbitrarily large number of particles can 
be stored in the condensate, no connection can be es- 
tablished between the number of particles and the heat 
capacity of the coolant. Thus condition (|13H cannot be 
considered as a constraint on Cb,o/Cf^ |2Cj . 

Recall that ST/T <C 1 was required to replace the finite 
differences in (JSJ) by derivatives. When the heat capac- 
ity of the Bose gas is initially much smaller that that 



of the Fermi gas, one finds that this condition is satis- 
fied for any value of rj, including rj = + where /(0) is 
finite and strictly positive. In this case, Eq. (|ll|l and 
its consequences apply also when one evaporates all the 
non-condensed atoms at each cooling cycle. 



IV. COOLING BY A BOLTZMANN GAS 



A Boltzmann gas has a total energy 

E b (N b ,T) = (S + 3/2)N b k B T. 



Equation Q then reduces to 



ST 
5N~ b 



[/(?,) - (S + 3/2)]T 



2aN f 



(5 + 3/2)N b 



(15) 



(16) 



The explicit dependence of the energy on Nb introduces 
a negative contribution in the numerator. However, the 
total numerator remains positive, since one can show, by 
integration by parts, that the condition 



/(r?) > 5 + 3/2 

is satisfied for all rj > 0. 

One finds that, by making a change of variables 



2a T N 



f 



C 



5 + 3/2 Tp N h C b 



(17) 



(18) 



the equation becomes separable. This leads to the con- 
stant of motion 

K = -——ln\X-l-x\ + —^—]iix-ln(N b /Nf), (19) 
A — 1 A — 1 



where 



A= 5T3f2 >0 - 



(20) 



The asymptotic behavior of the temperature in the limit 
of the complete evaporation of the Boltzmann particles, 
i.e., Nb — ► 0, falls into two regimes. The first is given by 
xq > A — 1, where Xq is the initial value of a;. In this case, 
the temperature tends to a finite non-zero value due to 
an insufficient number of Boltzmann particles. The ratio 
of the final to initial temperature is 



Hi 



rial 



T 



l-(A-l) 



-lV(A-l) 



Cf,o 



(21) 



The second regime is given by xq < A — 1, which implies 
A > 1. The inequality x<j < A — 1 is equivalent to the 
requirement that there be a minimal initial number of 
Boltzmann particles: 



c f,o 



k B [f(tj)-2{5 + 3/2)) 



(22) 
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FIG. 3: Shown is the temperature of a degenerate Fermi gas 
subject to sympathetic cooling by a Boltzmann g cLS ELS ct func- 
tion of the number of Boltzmann particles. Three cases are de- 
picted: Cf.o/Cb,o = (A — l)/2 (solid curve); C/^/C^o = A — 1 
(dashed curve); and Cfp/Cbfi = 2 (A — 1) (dot-dashed curve). 
In the latter case the number of Boltzmann particles is insuf- 
ficient to cool to T = 0. Here the choice of parameters was 
?7 = 6 for the energy cut-off parameter, Tq/Tf — 0.2 for the 
initial temperature, and a harmonic trapping potential. 



Note that A > 1 implies that the denominator is positive. 
Equation may be equally expressed as a requirement 
on the initial heat capacities of the fermions and Boltz- 
mann particles: 

C f , <C M (A-1). (23) 
In this case the temperature tends to zero as 

Tf \N f J { ' 

In Fig.[3]are shown three sample trajectories of Eq. (|16f> . 

For xq > A — 1 the final temperature is non-zero when 
all of the Boltzmann particles have been evaporated. For 
xq < A — 1 the final temperature is zero and all of the 
Boltzmann particles are used in the evaporation. A third 
curve shows the case xq = A — 1, which divides the two 
regimes: this corresponds to constant x, which implies 
T/T F oc N b /N f . 



Note that, for a harmonic potential, f{rj) can be calcu- 
lated exactly for a Boltzmann gas 01 ■ The critical value 
of if which corresponds to A — 1 = is then r\ c = 4.59 . . .. 
For rj < rj c the final temperature of the fermions after 
evaporation of all the Boltzmann particles is always non- 
zero. 



V. CONCLUSION 

We have demonstrated that, in an idealized experiment 
with no particle losses besides those due to evaporation 
of the coolant, sympathetic cooling of degenerate Fermi 
gases is not limited by the heat capacity of the coolant 
when the latter is a Bose gas with a condensate present. 
In fact, it is advantageous to obtain Cf/Cb 3> 1 in order 
to minimize the number of bosons that must be evapo- 
rated, where CV & are the heat capacities of the fermions 
and bosons. In this case, T = may be obtained provided 
that the initial number of bosons is sufficiently large: 

N bt0 > 2.7 N f T /T F , (25) 

where Nf is the number of fermions, To is the initial tem- 
perature which is much less than the Fermi temperature 
TV, and 77 = 6, a typical value of the energy cut-off for 
evaporative cooling, was chosen. 

In contrast, in the case of a Boltzmann coolant, the 
final temperature of the fermions depends on the ratio of 
the initial values of the heat capacities Cf^/Cb^- For 
typical values of the evaporation cut-off parameter 77, 
C/,o/C&,o must be smaller than a critical value on the 
order of unity to reach T = after evaporation of all the 
Boltzmann particles. 

Acknowledgments: We would like to thank Thomas 
Bourdel and Jean Dalibard for useful discussions, and 
Zoran Hadzibabic and Roland Combescot for useful com- 
ments on the manuscript. This work was supported by 
NSF grant no. mps-drf 0104447 Laboratoirc Kastler 
Brosscl is a research unit of l'Ecolc Normalc Supcricurc 
and of l'Universite Pierre et Marie Curie, associated with 
CNRS. We acknowledge financial support from Region He 
de France. 



[*] Present address: JILA, National Institute of Standards 

and Technology and Physics Department, University of 

Colorado, Boulder, CO 80309-0440 
[1] A. G. Truscott, K. E. Strecker, W. I. McAlexander, 

G. Partridge, and R. G. Hulet, Science 291, 2570 (2001). 
[2] F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, 

T. Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. 

Lett. 87, 080403 (2001). 
[3] G. Roati, F. Riboli, G. Modugno, and M. Inguscio, Phys. 

Rev. Lett. 89, 150403 (2002). 



[4] Z. Hadzibabic, C. A. Stan, K. Dieckmann, S. Gupta, 

M. W. Zwierlein, A. Gorlitz, and W. Ketterle, Phys. Rev. 

Lett. 88, 160401 (2002). 
[5] Z. Hadzibabic, S. Gupta, C. A. Stan, C. H. Schunck, 

M.W. Zwierlein, K. Dieckmann, and W. Ketterle, Phys. 

Rev. Lett. 91, 160401 (2003). 
[6] R. Onofrio and C. Presilla, Phys. Rev. Lett. 89, 100401 

(2002) . 

[7] C. Presilla and R. Onofrio, Phys. Rev. Lett. 90, 030404 

(2003) . 



5 



[8] K. Huang, Statistical Mechanics, 2nd ed. (John Wiley & 

Sons, New York, NY, 1987). 
[9] It is assumed for simplicity that the heat capacities are 
constant in the range of temperatures considered. 
[10] H. F. Hess, Phys. Rev. B 34 3476 (1986). 
[11] O. J. Luiten, M. W. Reynolds, and J. T. M. Walraven, 

Phys. Rev. A 53 381 (1996). 
[12] K.B. Davis, M.O. Mewes, W. Ketterle, Appl. Phys. B 60, 
155 (1995). 

[13] M. Wouters, J. Tempere, and J. T. Devreese, Phys. Rev. 
A 66 043414 (2002). 

[14] A gas of identical particles is in fact described by either 
Bose or Fermi statistics as the temperature approaches 
zero. However, for the sake of highlighting the effect of 
Bose statistics in contrast to classical statistics in the 
coolant, and in order to explain the meaning of the classi- 
cal thermodynamic argument mentioned above, we shall 
imagine that Boltzmann statistics can be applied down 
to zero temperature. 

[15] We do not claim here that zero temperature can be 
achieved in a real experiment: it would require an infinite 
time, as the collision rate between the fermions and the 
coolant tends to zero for T — > 0. Furthermore the semi- 



classical approximation used in this paper fails when the 
temperature drops below the energy of the first excited 
state of a particle of the coolant. 

[16] C. Cohen- Tannoudji, Cours de physique atomique et 
moleculaire, http:/ /www.lkb.ens. fr/~cct| 1997. 

[17] Note that such a density of states corresponds [lrt to 
a trapping potential of form V{x,y,z) = a|a;| 1/ ' 5ie + 
b\y\ 1/6y + c\z\ 1/s - where 5 X + 5 y + 8 Z = 8. 

[18] In the case of harmonic trapping, 5 = 3/2 is half integer 
so that f(r/) has a rather simple expression. For a Bose 
gas below T c we find f{rf) = [7? 3 gi(e"") + 3n 2 ^(e - ") + 
6r,g 3 (e^) + 6g 4 (e^)]/[i 1 2 g 1 (e-v) + 2r,g 2 (e^) + 2g 3 (e-^] 
where the Bose function is g a {z) = ^2 k>1 z k /k a ■ For a 
Boltzmann gas f(n) = [rf + 3r/ 2 + 6r/ + 6] /fa 2 + 2rj + 2]. 

[19] Handbook of Mathematical Functions, edited by M. 
Abramowitz and I. A. Stegun (National Bureau of Stan- 
dards, Washington, D. C, 1964). 

[20] One may wonder if the condition 11411 automatically im- 
plies (11311 . It turns out to be the case only if 77 is large 
enough, that is if (5 + 5/2)C/, < {[f(v)/f(0)] - [{S + 
5/2)/((5 + 3/2)]}C M . 



